Fidelity susceptibility for SU(2)-invariant states 
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We study the fidelity susceptibility of two SU(2)-invariant reduced density matrices. Due to 
the commuting property of these matrices, analytical results for reduced fidelity susceptibility are 
obtained and can be applied to study quantum phase transitions in SU(2)-invariant systems. As 
an example, we analyze the quantum criticality of the spin-one bilinear-biquadratic model via the 
fidelity approach. 
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Introduction — In 1984, Peres introduced the concept 
of fidelity to characterize quantum system responses to 
a perturbation [1]. It is of fundamental importance 
when studying quantum dynamics and has been applied 
to characterize two important phenomena in condensed 
matter theory, quantum chaos and quantum phase tran- 
sitions (QPTs) [1, i, i, i, i, 0, ij. It also becomes an 
useful concept in quantum information theory [j|, and 
has been used in study of quantum entanglement the- 
ory [l(]| , quantum teleportation [Til ] , and transformation 
of unknown states [12] etc. 

As an indicator of QPTs, various kinds of fidelity 
has been used in investigating quantum phase transition 
point, such as Loschmidt echo [H, ground-state fidelity 
[13|], the fidelity of first excited state 14 1, operator fi- 
delity [H Hi], and reduced fidelity Q][uM etc. The 
fidelity susceptibility (FS) ||, as the leading term of the 
fidelity, can be conveniently used to detect QPTs for its 
independence of the concrete values of small perturba- 
tions. 

Let us briefly introduce quantum fidelity and fidelity 
susceptibility. For pure states, fidelity is the absolute 
value of an overlap of two wave functions. One important 
case is the fidelity between the ground state \i/jq(x)) of 
Hamiltonian H(x) and a slightly different one \ifj Q (x + 5)), 

F= | <V>o(aO|Vo(:c + *)>!, (1) 
where S is a small deviation. Substituting the expansion 
\yj(x + S)} = \i/>(x)) + 5\4ix)} + S 2 /2\^(x)} + 0(S 3 ) (2) 
into the above equation leads to 



F = 1 + T ««o) + (MM + 2|(^o)| 2 ) 
= l-^((V'olV'o)-|(^ol^}| 2 ). 



(3) 



One may further define the FS as [f| 
2(1 -F) 



s 2 

= ^KVv,|V>o} 



= <Vo|^o> - \(MM\' 



(4) 



So, the FS is explicitly written out in terms of eigenstates. 
However, for mixed-state case, the corresponding fidelity 
and FS is relatively difficult to be achieved. One can use 
Uhlmann's fidelity [l!| 



F = tiVe 1/2 ^ 1/2 



(5) 



for two mixed states g and g and the corresponding FS 
can also be defined as above. In what follows, we analyze 
fidelity of SU(2)-invariant mixed states by this definition. 

For a many-body quantum state, by tracing out other 
degree of freedom but two particles, we have a two- 
particle reduced-density matrix which is generally a 
mixed state. Fidelity between reduced density matri- 
ces is called reduced fidelity [13] ■ For some interesting 
physical models such as the spin-one bilinear-biquadratic 
model [20l |2ll [22l . I23I ] and spin- half frustrated model, the 
reduced-density matrix displays a SU(2) symmetry. In 
this paper, we will study the fidelity of SU(2)-invariant 
state. The symmetry in this state faciliates greatly our 
study of fidelity and analytical results are obtained for 
the fidelity susceptibility. We also give an application 
of the results to study quantum phase transition in the 
bilinear-biquadratic model. 

SU (2) -invariant states and FS — Before proceeding, 
we make it clear that if a multi-spin state p dis- 
plays a global SU(2) symmetry ([p, J] = 0), the two- 
spin reduced-density matrix also has a SU(2) symmetry 
([/O12J1 +J2] = 0), where J = ji + ... + jjv is the col- 
lecting spin operator. The proof is straightforward. The 
commutator [p, J] = means that 



-J2J 



[J3 



-jn,p\- 



(6) 



After tracing out degree of freedom of spins 3 — > N, we 
have 
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[Pl2,jl + h] = TT3^Jv|j3 + ... +jn,p] = 0. 



(J) 
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An SU(2)-invariant state of two spins ji and j 2 can be 
written in the general form 



h+32 

Eaj 
2JTT 

J=\3i-h\ 



Pj, 



(8) 



where aj > 0, J2j a J = 1) an( l i s tne projector of 
spin- J subspace. Obviously, the density operator has 
eigenvalues aj/(2J + 1) with degeneracy 2J+ 1. 

One key observation from the above equation is that 
two different SU(2)-invairant density matrices p and p 
commute with each other. Thus, they can be diagonal- 
ized simutaneously, and the fidelity between them are 
given by. 



F 



3132 

E 

fc=i 



(9) 



4A fc 



(11) 



In deriving the above equation, we have used J^. A,; = 1 
and J^j d a Xi — J2i ^qAj = 0. Therefore, accordingto the 
relation between fidelity and FS F = 1 — x<5 2 /2 the 
FS xf corresponding to the matrix p is obtained as 



Xf 



E 



(d x x k y 

4A fc 



(12) 



This expression of fidelity susceptibility is valid for any 
commuting density matrices. It depends on nonzero 
eigenvalues of p and their first-order derivatives. 

Applying the Eq. (fl"2"|) to the SU(2)-invariant state p 
©, one obtains the FS for p as 



Xf 



31+32 I r\ n9 

J =32 -h 



4a 



J 



(13) 



where we assumed j 2 > ji without loss of generality. 
Now, we consider the following case of ji = 1/2 and 
j 2 > 1/2. As a j2 -i/2 + "33+1/2 = 1, Eq. flTSJ) reduces to 



XF 



(d x a h _ 1/2 ) 2 



4a j2 _ 1/2 (l-a j2 _ 1/2 ) 



(14) 



Parameter a j2 _i/ 2 can be written in terms of expectation 
of Heisenberg interaction on p, (ji • j 2 ), i.e., [2J] 

a j 2 -i/2 = 2 ^ + 1 (h - 2(ji • j 2 ))- (15) 
Thus, Eq. (fT4|) can be reexpressed in the following form 



XF 



0*2 - 2(ji - J2»C72 + 1 + 2(ji - ja» 



(16) 



We see that for the SU(2)-invariant state, the FS is com- 
pletely determined by the expectation value of Heisen- 
berg interaction and its first-order derivative. If we con- 
sider the case of two qubits, then the above equation 
reduces to 



XF 



4(9 K (ji-J2}) 2 



(17) 



where A^'s and A^'s are the eigenvalues of p and g, re- 
spectively. Since zero eigenvalues have no contribution 
to F, we only need to consider the nonzero ones. In the 
following, the subscript k in only denotes nonzero 
eigenvalues of p. 

Now we calculate fidelity of two slightly different den- 
sity matrices p(x) and p(x + 6) as a function of parameter 
x, where 5 is a small change of x. It is noticed that, for 
a small change S, we have 

\ k (x + 6)~\ k + (d x \ k )6+(d 2 x \ k ) 6 2 /2 + (5 3 ). (10) where 

Substituting this expression into Eq. © leads to the fi- 
delity given by 



(l-4(ji-j 2 ))(3 + 4(j 1 .j 2 ))' 

1 via a different 



which is just the FS obtained in Ref. 
approach. 

Now, we study the case of two qutrits, i.e., two spin 
ones. From Eq. ([8|), one has 

a = (Pa) = (P12), 



ai = (P L ) = ~(l-(S i3 )), 



a 2 = (P 2 ) = -(l-2(P 12 ) 



<Si2», 



(18) 



P12 = gCOi • J2) 2 - 1], 
5*12 = ji h + (ji ■ h) 2 - 1, 



(19) 



are singlet projection operator and swap operator, re- 
spectively. Substituting Eqs. (TTg|) and (Til?)) into Eq. |T 
leads to the FS for two qutrits 



1 

XF = ~ 4 



(d x (p 12 )f , (d x (s 12 )f 



(P 12 ) 2(1 - (S l2 )) 

{d x (s 12 )-2d x (p 12 )f- 



2(1 - 2(P 12 ) + (5i2» 



(20) 



The FS is detrmined by two expectation values (Pl 2 ) and 
(S'12) and their first-order derivatives. Below, we will ap- 
ply this formula to the study of the bilinear-biquadratic 
model. 

Applications to spin-one systems — Spin Heisenberg 
chains attract more attention since Haldane predicted 
that the one-dimensional chain has a spin gap for inte- 
ger spins [1^1 . In these studies, the bilinear-biquadratic 
model has played an important role [2(J, HH H2, HH . The 
corresponding Hamiltonian is given by 



TV 



(21) 



Hbb = ^cos6»(j 4 -jj+i) +sin6»(j l ■ j J+ i) 2 , 

fcos 9Si,i+i + 3(sin 9 - cos 0)P i)i+1 ] + N sin 9 



i=l 
N 

Ei 

i=l 
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In deriving the last equality, we have used Eq. (|T9]) . Here, 
ji denotes spin-1 operator at site i, and we have assumed 
the periodic boundary conditions. The Hamiltonian ex- 
hibits a SU(2) symmetry, and displays very rich quantum 
phase diagrams (27| . 

From Hellmann-Feymann theorem for ground state, 
one can easily find that 

(P12) = i(sin6>e + cos6»e - 1), 
(S12) 

— (cos 9 + sin 0)eo + (cos 9 — sin 0)e o — 1, (22) 

and the their first-order derivatives 
cos9 

(-P12) = -g-(eo + e ), 
(Sx 2 )'= (cos0-sin0)(e o + e o '), 



(s 12 y - 2{p 12 y = 



cos 9 



-Bin0) (eo + eft). (23) 



Here, eo denotes the ground-state energy per site. Sub- 
stituting the above two equations into ([2U)l . one obtains 
the FS in terms of eo, e , and e ' as follows 



XF 



(eo 



4f 



cos 2 9 



3(sin0eo + cos6>e — 1) 



+ 



2 [2 - (cos6» + sin0)e o - (cosfl - sin6»)e ] 
(cos6> - 3sin6>) 2 
+ 6 [2 + (3 cos 8 + sin 9)e Q + (cos 9-3 sin 9)e' ] 

(24) 

One key observation is that the numerators of the 
above two expressions happen to be proportional to 
(eo + e ') 2 . Then, we infer that if the second deriva- 
tive of the ground-state energy is singular at the critical 
point, the FS is singular too. On the other hand, it is 
known that the divergence of the second derivative of 
the ground-state energy reflects the second-order QPTs 
of the system, which is shown in Ref. [l3[ explicitly as 



dle 



j^ 2\(V n \d a H\^ r , 



N(E - E n ) 



where \^ n ) is the eigenvector corresponding to the eigen- 
value E n . It shows that the vanishing energy gap in the 
thermodynamic limit can lead to the singularity of the 
the second derivative of the ground-state energy. There- 
fore, the two-spin FSs can exactly reflects the second- 
order QPTs of the global system in this model. 

We use the exact-diagonalization method to calcuate 
the ground-state energy and then numerical results of FS 
is obtained from Eq. lfM]) . In Fig. 1, we plot the FS as 
a function of 9 for a sample of 12 spins. We observe a 





0.5 1 
8/71 



FIG. 1: Ground-state fidelity susceptibility as a function of 
9/tt in the bilinear-biquadratic model. 



occurs. For 7r/2 < 9 < 5ir/4, the ground state is ferro- 
magnetic and degenerate. In this range, the FS is zero. 
However, at 9 = 7ir/4, corresponding to a QPT point 
separating dimerized phase (5ir/4 < 9 < 7tt/4) and Hal- 
dane phase, one cannot find any anomalous behaviors of 
the FS. 

Conclusions — We have studied the FS in SU(2)- 
invariant states and have obtained exact analytical ex- 
pression of the FS for any spins j\ and j 2 . This implies 
that the results are applicable to not only equal-spin but 
also mixed-spin systems. Furthermore, one can use the 
FS to study properties of SU(2)-invariant physical sys- 
tems in a finitc-tcmpcrature thermal state. As an ap- 
plication, we have studied relations between the FS and 
QPTs in the bilinear-biquadratic model. For this model, 
one can infer that the two-spin FS can exactly reflects 
the second-order QPTs of the global system. Here, we 
restrict us to study SU (2) -invariant states of two spins, 
after tracing out other spins of a many-body states. One 
challenge for further investigation is to study ./V-spin 
(N > 3) SU(2)-invariant reduced density matrix. 
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